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Abstract This contribution reviews a selection of findings on atomic density func- 
tions and discusses ways for reading chemical information from them. First an ex- 
pression for the density function for atoms in the multi -configuration Hartree-Fock 
scheme is established. The spherical harmonic content of the density function and 
ways to restore the spherical symmetry in a general open-shell case are treated. The 
evaluation of the density function is illustrated in a few examples. In the second part 
of the paper, atomic density functions are analyzed using quantum similarity mea- 
sures. The comparison of atomic density functions is shown to be useful to obtain 
physical and chemical information. Finally, concepts from information theory are 
introduced and adopted for the comparison of density functions. In particular, based 
on the Kullback-Leibler form, a functional is constructed that reveals the periodic- 
ity in Mendeleev's table. Finally a quantum similarity measure is constructed, based 
on the integrand of the Kullback-Leibler expression and the periodicity is regained 
in a different way. 



1 Introduction 

Density Functional Theory (DFT) plays a prominent role in present day investiga- 
tions of the electronic structure of atoms and molecules. Within DFT the electron 
density function plays a central role as it caries all the information to describe the in- 
vestigated system. The idea that all physical and chemical information is contained 
in the density incited the present authors to try and recover some of this informa- 
tion. Although the proof of the Hohenberg-Kohn theorems 1 1 1, which guarantee the 
presence of all physical information in the density function, is generally said to be 
disarmingly simple, it does not provide a method to get to the relevant information. 
The continued search for improved energy functionals is the most evident example 
to illustrate the challenge researchers are confronted with. 

In recent years, the present authors have developed an interest in obtaining chem- 
ical information from atonic density functions. The application of concepts from 
quantum chemistry shows that some particular aspects of physical and chemical in- 
terest can be read from the density functions. In particular the comparison of density 
functions using quantum similarity measures or functionals from information theory 
plays an important role. The original goal of the work was to find a way of regaining 
the periodicity in Mendeleev's table through the comparison of density functions. 

The purpose of this contribution is to give an overview of the results which center 
around the atomic density function and the recovery of the periodicity. Since all the 
calculations are based on atomic density functions, it is appropriate to revisit the 
construction of these densities in some depth. First a workable definition of the 
density function is established in the framework of the multi-configuration Hartree- 
Fock method (MCHF) 1 2 1 and the spherical harmonic content of the density function 
is discussed. A spherical density function is established in a natural way, by using 
spherical tensor operators. The proposed expression can be evaluated for any multi- 
configuration state function corresponding to an atom in a particular well-defined 
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State and a recently developed extension of the MCHF code IS is used for that 
purpose. Three illustrative examples are given. In the next section relativistic density 
functions for the relativistic Dirac-Hartree-Fock method [4| are defined. The latter 
will be used for a thorough analysis of the influence of relativistic effects on electron 
density functions later on in this paper 

The analysis of atomic density functions can be furthered by comparing them in 
pairs. Specifically, the use of quantum similarity measures and indices as defined 
by Carbo Q has shown that particular influences on the density functions can be 
estimated in this way. Here this feature is demonstrated by reviewing three case 
studies: i) the L5-term dependence of Hartree-Fock densities, ii) the comparison of 
atoms throughout the periodic table [6| and iii) the quantitative evaluation of the 
influence of relativistic effects, via a comparison of non-relativistic Hartree-Fock 
densities with Dirac-Hartree-Fock relativistic densities Q. 

In the final part of this contribution, information theory is introduced. After a 
brief revision of the relevant concepts, a functional based on the KuUback-Leibler 
measure [8 1 is constructed for the investigation of atomic density functions through- 
out Mendeleev's table. Since the quantum similarity does not reveal the expected 
periodic patterns, it is significant to show that it is actually possible to regain the pe- 
riodicity by constructing an appropriate functional |6|. By considering the integrand 
of the Kullback-Leibler measure and comparing it locally for two atoms, a quantum 
similarity measure can be constructed which does reveal periodic patterns Q. 



2 The multi-configuration many-electron wave function 

The total energy which results from the Hartree-Fock equations is due to electrons 
moving independently in an averaged, central potential. Any improvements to the 
energy (still in the context of the non-relativistic Schrodinger equation) are said 
to be due to correlation effects as a direct consequence of the electron-electron 
interaction. In fact it is common to define the correlation energy as 

Eqorr = -Eexact - EuF , (1) 

where £exact is the exact non-relativistic energy and EfjF the energy due to the 
solutions of the corresponding Hartree-Fock equations. 

A particularly natural way to improve the Hartree-Fock energy - i.e. include cor- 
relation energy - is by departing from the single-configuration approximation. In 
Hartree-Fock calculations a rigid orbital picture is assumed, where electrons have a 
fixed place in a given electron configuration. By allowing electrons to occupy dif- 
ferent orbitals and allowing several electron configurations, the variational approach 
can be applied on a significantly larger set of trial wave functions. 

In the multi-configuration Hartree-Fock (MCHF) approach, the A^-electron wave 
function ^aLSM^Ms is a linear combination of M configuration state functions (CSFs) 
^aiLSM^Ms which are eigenfunctions of the total angular momentum L^, the spin 
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momentum and their projections and S^, with eigenvalues h^L{L+ 1), h^S(S + 
1), fiM]^ and fiMs, respectively 

M 

•fttLWiMslxi, • --xn) = <PiaiLSMLMs;xi,- --xn) , (2) 

1=1 

with 

M 

(3) 

/=1 

The mixing coefficients {c,} and the radial functions {/?„./. (r)}, constituting the 
one-electron basis, are solutions of the multi-configuration Hartree-Fock method 
in the MCHF approach. For a given set of orbitals, the mixing coefficients may 
also be the solution of the configuration interaction (CI) problem. The relativistic 
corrections can be taken into account by diagonalizing the Breit-Pauli Hamiltonian 
||9l in the L57-coupled CSF basis to get the intermediate coupling eigenvectors 

M 

1'aJMixi , • • • xa,) = 2^ a/ <P{aiLiSiJM; xi,---xn) , (4) 

1=1 

with 

M 

Lkl'^l. (5) 

1=1 

In the MCHF approach, the trial wave functions are of the form Q and the energy 
expression for a given state becomes 

M M 

'■ j 

where 

ah abcd;k 

where a = [nJa) and b = [iihlb) and the sums run over all occupied orbitals in the 
respective configuration / and j. 

In the MCHF context the variational principle is applied to the energy functional 
in equation (|6|. A stationary solution is obtained by minimizing the energy with 
respect to variations in the radial wave functions Pni{r) = rR„i{r) satisfying the 
orthonormality conditions 

Nnl.n'l = [ Pnl{r)Pn'l{r) dr ^ 5„ni (8) 

Jo 

and 

M 

Yc]^l. (9) 

1=1 
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The energy expression (|6]l can be written in matrix notation 

E{aLS) = C'HC , 



(10) 



where H = {Hjj), C is the column matrix of the expansion coefficients and C its 
transpose. 

This gives the eigenvalue problem for the expansion coefficients 



To obtain a self consistent field solution to the MCHF problem, two optimiza- 
tions need to be performed i.e. one for the variation of the one-electron radial or- 
bitals {P„;(r)} in the wave function and one for the expansion coefficients. This 
can be done by consecutively iterating, first the orbital optimization followed by the 
coefficient optimization. 



3 On the symmetry of the density function 

In the first part of this paragraph we review the study of the spherical harmonic 
content of the density function for atoms in a well defined state (|2]i or Q. The 
spherical density functions, which reveal the familiar shell structure, are discussed 
and illustrative examples are given. 



3.1 The non-spherical density function 

The so-called "generalized density function" ifTOl or the "first order reduced density 
matrix" ifTTIl is a special case of the reduced density matrix lfT2l[T0ll 



7l(xi,x'l) =A' / f (xi,X2,...,XAr) "f^* (x' i , X2 , . . . , Xa,) d^2---d^N , (12) 



where 'f'(xi ,X2, . . . ,XAr) is the total wave function of an electron system and 
f *(xi,X2, . . . ,XAf) is its complex conjugate. The spin-less total electron density 
function p(r) is defined as the first order reduced density matrix, integrated over 
the spin and evaluated for xi = x'l 



This electron density function is normalized to the number of electrons of the system 



HC = C£ . 



(11) 




(13) 




(14) 
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As discussed in ifTTI . the single particle density function can be calculated by eval- 
uating the expectation value of the 5 (r) operator, 

p(r) = J ^{xi,X2,...,xn) d{r)^*{xuX2,...,XN) dxidx2...dxN , (15) 

where 5 (r) probes the presence of electrons at a particular point in space and can 
be written as the one-electron first-quantization operator 

5(r) = f 5(r-r,). (16) 

!=1 



The exact spin-less total electron density function ( 15 1 evaluated for an eigenstate 
with well-defined quantum numbers {LSM^Ms) 



1 ^ 



(17) 



Im ^ i= 1 



becomes 



where 



p(^)«M,M, ^ ^p(,)«^^.M.y2/o(T?,(p) , (18) 

/=o 



" r2 ^ ' \~Ml Ml 

N 

X (faisll Yi,{&i,(Pi)\\1'aLs) ■ (19) 

i=l 

This result, which can be found by applying the Wigner-Eckart theorem llT3l . recov- 
ers Fertig and Kohn's analysis [14] for the density corresponding to a well-defined 
[LSMiMs] eigenstate of the Schrodinger equation. The spherical harmonic compo- 
nents in the density are limited to /-even contributions, because the bra and the ket 
need to be of the same parity % = (—1)^' '' . 

In this paper, the authors observed that the self-consistent field densities obtained 
via the Hartree and Hartree-Fock methods generally violate the specific finite spher- 
ical harmonic content of p (r)^'^'^'^'^-^ . They also mention that this exact form can be 
obtained by spherically averaging the effective potential, yielding single-particle 
states with good angular momentum quantum numbers. The atomic structure soft- 
ware package ATSP2K |2| applies this approach, as was done in the original atomic 
Hartree-Fock theory 1,15.,16..17 |. This implies two things: i) the density function 
p ^^y^SMlMs calculated from any multiconfiguration wave function of the form (j2]), 
is not a priori spherically symmetric, ii) this density function will contain all spher- 
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ical harmonic components (up to 2L) as long as the one-electron orbital active set 
spanning the configuration space is /-rich enough. 

The density fiinction can also be expressed in second quantization |10|. Intro- 



ducing the notation q = n^lqini^ms^^ for spin-orbitals, expression ( 12 1 becomes 



7i(xi,x'i) i/;(x'i)v/^(xi) , (20) 



where Dpq are elements of the density matrix which are given by 

Dp„ = {'¥\ala^\'¥) . (21) 
The sum in eq. (|20]l runs over all possible pairs of quartets of quantum numbers 



p and q. The spin-less density function ( 13 1 calculated from p(r) = (f |5(r)|f 



using the second quantized form of the operator 



yields 



p(r) -^D,, 5,„,^,,„,^^ 7?,;,/r)F,; (tJ,(p)7?„,^,(r)r,^,„,^(t?,(p) . (23) 



To illustrate the spherical harmonics content of the density in the Hartree- 
Fock approximation, consider the atomic term \s^2p^{ ^P)3d for which the 
{Ml, Ms) = (+3, +3/2) subspace reduces to a single Slater determinant 



faLSM.M, = 'I>{ls^2p^{ ^P)3d V+3,+3/2) = \lsls2p+i2po3d+2\ . (24) 



When evaluating ( |23| , all non-zero D^^-values appear on the diagonal (p — q), yield- 
ing 

p (r)'^«.+3/2 ^ I v/i,(r) |2 + 1 v/Tj(r) |2 + 1 y/jp^, (r) p + 1 y/j^^ (r) ^ + 1 {r)\\ (25) 

This density has a clear non-spherical angular dependence. However, referring 
to US) 

J J 

w}j^{^) = \YjM{-&,(p)\^^Y.^„{JM)P2n{cQS-&)^Y.^'n{JM)Y2„0{-&M (26) 

n=0 n=0 

one recovers the even Legendre polynomial content of the density, although not 



reaching the (2L = 6) limit ygoC??, <p) of the exact density ( 18 1. This limit will be 
attained when extending the one-electron orbital active set to higher angular mo- 
mentum values for building a correlated wave function. 
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As discussed in detail in fS), contributions to the density function corresponding 
to {p ^ q) can appear through off-diagonal matrix elements in the CSF basis. These 
contributions will be present for parity conserving single electron excitations of the 
type \hq) Ihq). 

The "offending" spherical harmonic contributions described by Fertig and Kohn 
fl4\ do not occur in the MCHF calculation of the density function, whatever the 
maximum /-value of the orbital active space |i3j|. 



3.2 The spherical density function 

A spherically symmetric density function can be obtained for an arbitrary CSF 
'PaLSMtMs by averaging the (2L + 1 ) {2S + 1 ) magnetic components of the spin-less 
density function 



where p (r)^'^*''-'*^^ is constructed according to eq. ( 23 i 



p(r)LW.M, =Y^{cp^^sM,Ms\ala,\<PaLSM,Ms) K,„n.., ¥;{r)W,ir) . (28) 
PI 



Applying equations (27 1 and (28 i for the atomic term Is 2p { P)3d F consid- 



ered in the previous section, we simply get 

P (r)'^ = ^ { 2Pl (r) + 2Pl (r) + pI (r) } , (29) 
which is, in contrast to eq. ([25|), obviously spherically symmetric. The sum over 



[Ml, Ms) performed in (27 i guarantees, for any «Z-subshell, the presence of all nec- 
essary components {y/„,, | mi — ~l, ... + 1} with the same weight factor, which per- 
mits the application of Unsold's theorem |fT9l 

|F„„,(i?,(p)|2 = ^ (30) 

m; =— / 

and yields the spherical symmetry. This result is valid for any single CSF 

p«" = iE?«'^3w, (31) 

nl 

where q„i is the occupation number of n/-subshell. Its sphericity explicitly appears 



by rewriting ( 3 1 1 as 
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p(r)-p(r) \Yooi&,(p)\' 



D{r) 



(32) 



with 



(33) 



and 



D{r) ^ r^p{r) ^ Y.^„iP^M = '"XV'-) 



(34) 



The radial distribution function D{r) represents the probability of finding an elec- 
tron between the distances r and r + dr from the nucleus, regardless of direction. 
This radial density function reveals the atomic shell structure when plotted as func- 
tion of r. Its integration over r gives the total number of electrons of the system 



D{r) dr 



r^p{r) dr = Y,qnl 

nl 



(35) 



Where above the spherical symmetry of the average density (27 1 is demonstrated 
for a single CSF thanks to Unsold's theorem, it can be demonstrated in the general 
case by combining (27 1, (18 i and the 3- j sum rule 1 13 1 



Ml 



L k L 

-Ml Ml 



(36) 



for each k = 21 contribution ( 19 1. However, the radial density p(r) will be more 
complicated than (33i, involving mixed contributions of the type P„//(r)P„;(r) = 
r^R„ii{r)R„i{r), as developed below. 

Instead of obtaining a spherically symmetric density function by averaging the 
magnetic components p (r)^™^'^^ through eq. ( 27 1, one can build a radial density 



operator associated to the function ( [34| which is spin- and angular-independent, i.e. 
independent of the spin (a) and angular {■&, (p) variables. Adopting the methodology 
used by Helgaker et al fTP\ for defining the spin-less density operator, we write a 
general first quantization spin-free radial operator 



in second quantization as 



i=l 



PI 



where /„„ is the one-electron integral 



(37) 



(38) 



fpq = / W*p{^)f{r)Wq{^)r^vii-&drd-dd^da . 



(39) 



Applying this formalism to the radial density operator 



12 Contents 

5(r) = f5(r-n), (40) 

1=1 

and using the spin-orbital factorization for both p and q quartets, we obtain the 
second quantization form 

Hr) =Y^dpc,{r)alaq , (41) 

with 

dp,,{r) = 8,ni^mi^ 8m,^ms^ KplpidRn.l.irV , (42) 

where the Kronecker delta arises from the orthonormality property of the spherical 
harmonics and spin functions. With real radial one-electron functions, the operator 
( |4T] i becomes 

8{r)= Yj Sl'l Sm'm, Sm',m,,al,,,^,^^,a„i,nim,,R„'i'{r)Rnl{r)r^ (43) 

I Jl I I I ^ ^ 

n ,/ ,mi,m\,n.L,mi ,m_^, 

= E E '^l'lm,m,'^nlm,m,Rn'l{r)Rnl{r)r^ ■ (44) 

Its expectation value provides the radial density function D(r) = r^p (r) = 4-71 r^p (r) 
defined by Q and p?] ). 

Building the coupled tensor of ranks (00) from the [2(2/ + 1)] components of the 
creation and annihilation operators ll20l 



/ t N (00) _ 1 ^ 



the operator ( 43 i becomes 



5{r) - -^v/2(2/ + l)i: (al./a„,)|^™^ 7;„,,(r)7;„,(r)r2 . (46) 

/ n' .n 

The expectation value of this operator provides the spherical density function for 



any atomic state. Note that, in contrast to (28 1, the tensorial ranks (00) garantee the 



diagonal character in L, S,Ml and Ms, thanks to Wigner-Eckart theorem 

{aLSMLMs\T^^^^\a'L'S'M'i^M's) - (-1)L+s-m,-Ms 

L L' \ f S OS' 



X (aL^llrC"" II a'L'5') . (47) 
Moreover, the Mi/Mg independence emerges from the special 3 j-symbol 
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'ifij rrij J 

In other words, where the non-spherical components are washed out by the aver- 
aging process ( 27 1, they simply do not exist and will never appear for the density 
calculated from (46 1, for any {Mi^Mg) magnetic component. 

The radial distribution function D[r) = r^p (r) can be calculated from the expec- 
tation value of the operator (|46]), using the wave function (j2]l or Q. In the most 
general case (expansion Q), using the (L5)7-coupled form of the excitation opera- 
tor, 

/ ^mO / x(00) 

(al/,a„/j^ = (^a„,,a„,J^^ , (49) 



one obtains 



with 



and 



{^aJM\d{r)\^aJM) = (-l/"^ ( _M m) ("^-jWFP'W^aj) (50) 



= -17^(^1 (a:,a„,)^''"''' /p («'/,«/) , (51) 

/—I n.n' 



Ip{n'lM){r) = R„.,{r)R„,{ry . (52) 

The diagonal reduced matrix element (RME) evaluated with the Breit-Pauli eigen- 
vector Q has the following form 



where the RME in the {LS)J coupled basis reduces to 

{0{aiLiSiJM)\\Fl''^^\^{ajLjSjJM)) 



r(00)| 



and 



X {^{aiLiSi)\\Fp\\0{ajLjSj)) (54) 

^3 = -Ev/2(2/+i)^ (««'/««')ir ("'^'"'^ • ^^^^ 

/=! n.n' 



From a comparison of the operator (55 i with the non-relativistic one-body Hamil- 
tonian operator (see eq. (A5) of 121]), one observes that the angular coefficients of 
the radial functions Ip {n'l,nl) (r) are identical to those of the one-electron Hamilto- 
nian radial integrals /„// „/ , as anticipated from McWeeny analysis 1 10|. These angu- 
lar coefficients can be derived by working out the matrix elements of a one-particle 
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scalar operator Fp between configuration state functions with u open shells, as 
explicitly derived by Gaigalas et al 1221 . 



4 Three tangible examples 

First the evaluation of the density function and the influence of correlation effects 
is illustrated by plotting in figure [l] the radial density distribution D{r) — r^p{r) 
for a CAS-MCHF wave function of the beryllium ground state (Be \s^2s^ ^S), us- 
ing a « = 9 orbital active set. In the same figure, the Hartree-Fock radial density 
is compared with the one obtained with two correlation models: i) the n — 2 CAS- 
MCHF expansion, largely dominated by the near-degeneracy mixing associated to 
the Layzer complex ls^{2s^ ^-Ip^) and ii) the n = 9 CAS-MCHF. From the plotted 
results we notice that the density of the n = 2 calculation already contains the ma- 
jor correlation effects, compared to the n = 9 calculation. Indeed, the density does 
not seem to change a lot by going from the n = 2 to the n = 9 orbital basis, the 
valence double excitation Is^lp^ contributing for 9.7% of the wave function. From 
the energy point of view however, this observation is somewhat surprising (see Ta- 
ble [T}: the correlation energy associated to the n = 2 CAS-MCHF solution "only" 
represents 47% of the n — 9 correlation energy. 




Fig. 1 Density of ' S Be ground state for different CAS-MCHF wave function as compared to 
Hartree-Fock (HF). Density differences have been scaled by a factor 100. 



As a second example, we illustrate the influence of relativistic effects - in the 
Breit-Pauli approximation - on the density function of the Be-like 0"*+ atom, by 
comparing the densities of the fine-structure states \s^2s2p ^Pq, ^P^ and ^P2- From 
the plots in Figure|2]and the data given in Table|2]we observe that the largest energy 
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model 


energy (a.u.) 


correlation energy (a.u.) 


HF 

;j = 2-CAS 
n = 9-CAS 


-14.573 023 
-14.616 856 
-14.667 013 


£?=2 = £«=2 - £«f = 0.043 832 
E'^=^^ = - E"^ = 0.093 986 



Table 1 Total energy for the ground state of Be with different correlation models. 



difference corresponds to the largest difference in density function. More bound is 
the level, higher is the electron density in the inner region. The influence of rela- 
tivistic effects on the density function will be discussed thoroughly below. 



model 


energy (a.u.) 


energy difference (a.u.) 


'2 


-68.032 086 
-68.031 473 
-68.030 102 


E{^Pi- ^Po) =0.000 613 
E^Pi- Vi) =0.001 370 



Table 2 Total energy for 2s2p fine-structure levels of 0"*+ 

Finally a third example is given that is relevant when studying the electron affini- 
ties, as it is often interesting to investigate the differential correlation effects be- 
tween the negative ion and the neutral system |23|. Figure [3] displays the density 
functions D{r) of both the [Ne]3i^3/:''* ground state of neutral Sulphur (S) and 
the [Ne]3i^3/7^ ground state of the negative ion S^, evaluated with elaborate 
correlation models [23], together with their difference AD{r). The latter integrates 
to unity and reveals that the "extra" electron resides in the valence area. 



5 Relativistic density functions 

5.1 Relativistic multi-configuration wave functions 

In the relativistic scheme, the atomic wave function is, in the most general case, a 
combination of configuration state functions 

\nJMj) ^Y.'^ylvnJMj) , (56) 

V 

eigenfunction of the inversion operator /, the total angular momentum and its 
projection J^. v denotes all the necessary information for specifying the relativistic 
configuration unambiguously. The CSFs are built on the one-electron Dirac four- 
spinor 

^.(A-U ^'W^^'^^) ^ (57^ 
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where Xicli^) is a two-dimensional vector harmonic. It has the property that 
K\jfi{r) = K\ifi{r) where = j3((7 - L + l). The large {P(r)} and small {Q[r)} com- 
ponents are solutions of a set of coupled integro-differential equations ( [62| ) |24|. 
The mixing coefficients {cy} are obtained by diagonalizing the matrix of the no-pair 
Hamiltonian containing the magnetic and retardation terms ll25l . The two coupled 
variational problems are solved iteratively. For a complete discussion on relativistic 
atomic structure we refer to |26|. 

It is to be noted that the relativistic scheme rapidly becomes more complicated 
than the corresponding non-relativistic one. For example, if the ground term of Car- 
bon atom is described, in the non-relativistic one-configuration Hartree-Fock ap- 
proximation, by a single CSF | \s^2s^2p^ ^P), the relativistic equivalent implies the 
specification of the 7-value. For J — corresponding to the ground level of Carbon, 
the following two-configuration description becomes necessary 



|"lsW2/?2"(y = 0)) = ci\lshs^{2p*)^{J = 0)> +C2\lshshp^{J = 0)) , (58) 



impUcitly taking into account the relativistic mixing of the two L5-terms ( '5 and 
^P) arising from the 2p^ configuration and belonging to the / = subspace. p* and 



p in expression (58 i correspond to the y-values, j 
(k = —2), respectively. 



1/2 (fc = +1) and j = 3/2 



5.2 Multi-configuration Dirac-Hartree-Fock equations 

For the calculations of relativistic density functions we used a multi-configuration 
Dirac-Fock approach (MCDF), which can be thought of as a relativistic version 
of the MCHF method. The MCDF approach implemented in the MDF/GME pro- 
gram im |27]| calculates approximate solutions to the Dirac equation with the effec- 
tive Dirac-Breit Hamiltonian ETll 

//^«(ri,r2,...,M==£/i^(r,) + i:/!^(r„r,) , (59) 

(=1 i<i 

with 

/!^(r,)=ca-p + c2(j3-l)\^(r) (60) 

and 

R, s 1 aiOLi , COs(a);;r;;) — 1 

h\ri,rj) = '-^ cos(a),;r,;) + (a • V),. (a • V) . ^ • (61) 

The total Hamiltonian contains three types of contributions: the one-electron 
Dirac Hamiltonion, the Coulomb repulsion and the Breit interaction. These contri- 
butions, which appear in the energy expression, give rise to radial integrals, which 
need to be calculated for the two component wave function. We will simply state 
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the MCDF equations, which can be obtained by applying the variational principle to 
the energy expression, for variations in the expansion coefficients c,- in equation ( 56 1 
and both the large and the small components of the radial wave function. The coef- 
ficients c, can be determined from the diagonalization of a Hamiltonian matrix and 
the radial components can be optimized by solving the coupled integro-differential 
equations, here given for the orbital A 



(62) 

where the summation over B contains only the contributions for Ka — Kb, where Va 
is the sum of the nuclear and the direct Coulomb potentials and where Xp^ contains 
all the two electron integrals, except the instantaneous direct Coulomb repulsion. 
With the Lagrangian parameters Eab the orthonormality constraint 

/ {PAir)PB{r) + QA{r)QB{r)} dr ^ (63) 

is enforced. 



^ + a -l + aVA{r) 



5.3 Relativistic density functions 

For the purpose of quantifying the relativistic effects on the electron density func- 
tions, which are discussed later on in this contribution, we evaluate Dirac-Fock 
density functions, using a point nucleus approximation. In this section we describe 
how density functions can be obtained. 
By averaging the sublevel densities 

+ ^) Mj = -J 

the total electron density becomes spherical for any open-shell system, as found in 
the non-relativistic scheme (see section[J!2|) and can be calculated from 



UK 



where q^^ is the occupation number of the relativistic subshell (njc). Expression ( 65 i 
can be considered as the relativistic version of equation (27 1. 
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6 Analyzing atomic densities: concepts from quantum chemistry 



6.1 The shape function 

In the context of information theory (cf. section[8]) the shape function, defined as the 
density per particle 

a(r)^§^, (66) 
where is the number of electrons, given by 



N = Jpir)dr, (67) 

plays a role as carrier of information. In particular, the shape function is employed 
as a probability distribution which describes an atom or a molecule. 

The shape function first came to the scene of quantum chemistry in 1983 with 
the work of Parr and Bartolotti ll28l . Although the shape function had appeared 
before in another context, it is due to Parr and Bartolotti's work that the quantity 
owes its name. Just as for the density function, the shape function can be shown 
to determine the external potential v(r) and the number of electrons ll29l and 
so completely determines the system. The Kato cusp condition [^O"! leads to the 
nuclear charges and the relationship between the logarithmic derivative of the shape 
function and the ionization potential determines the number of electrons ||29l . On 
this basis a Wilson-like argument has been constructed 1311 (similar to the original 
DFT), confirming the shape function "as carrier of information" |32 1. Relationships 
between the shape function and concepts from conceptual DFT were established. 
In f33l a slightly different perspective is given on the fundamental nature of the 
shape function. 



6.2 Quantum similarity 

In chemistry the similarity of molecules plays a central role. Indeed, comparable 
molecules, usually molecules with a similar shape, are expected to show similar 
chemical properties and reactivity patterns. Specifically, there chemical behavior is 
expected to be similar t34j . The concept of functional groups is extensively used 
in organic chemistry 1351 . through which certain properties are transferable (to a 
certain extent) from one molecule to another, and the intense QSAR investigations 
in pharmaceutical chemistry 1 36 1 are illustrations of the attempts to master and ex- 
ploit similarity in structure, physicochemical properties and reactivity of molecular 
systems. 

Remarkably, the question of quantifying similarity within a quantum mechani- 
cal framework has been addressed relatively late, in the early 1980's. The pioneering 
work of Carbo and co-workers ||5]|32l led to a series of quantum similarity measures 
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(QSM) and indices (QSI). These were essentially based on the electron density dis- 
tribution of the two quantum objects (in casu molecules) to be compared. The link 
between similarity analysis and DFT |38 39 1 built on the electron density as the 
basic carrier of information, and pervading quantum chemical literature at that time, 
is striking. 

The last 15 years witnessed a multitude of studies on various aspects of quantum 
similarity of molecules (the use of different separation operators 1371 . the replace- 
ment of the density by more appropriate reactivity oriented functions 1381 [39l within 
the context of conceptual DFT |40|, the treatment of enantiomers l 3ni4Tll42ll43l ). 
With the exception of two papers by Carbo and co-workers, the study of iso- 
lated atoms remained surprisingly unexplored. In the first paper l44l atomic self- 
similarity was studied, whereas the second one f45l contains a relatively short study 
on atomic and nuclear similarity, leading to the conclusion that atoms bear the high- 
est resemblance to their neighbours in the Periodic Table. 

The work discussed below is situated in the context of a mathematically rigorous 
theory of quantum similarity measures (QSM) and quantum similarity indices (QSI) 
as developed by Carbo l37l l5l. Following Carbo, we define the similarity of two 
atoms (fl and b) as a QSM 

Z„i(i2) - J p„(r)i2(r,r')Pi(r') drdr' , (68) 

where £2 (ri , is a positive definite operator. RenormaUzation to 

sin = ' , (69) 

yields a QSI SIq with values comprised between and 1 . 

The two most successful choices for the separation operator i2(r,r') are the 
Dirac-delta 5(r,r') and the Coulomb repulsion jTzpy • Ths first is known to reflect 
comparison of geometrical shape of molecules, whereas the second is said to reflect 
the charge concentrations l44l . 



7 Analyzing atomic densities: some examples 

In the previous sections the calculation of the density function was discussed and a 
methodology for comparing them was introduced. In this section a quantitative anal- 
ysis of atomic density functions is made. It seemed interesting to employ concepts 



from molecular similarity studies (cf. section 6.2 1, as well developed in quantum 
chemistry and chemical reactivity studies. Here molecular quantum similarity mea- 
sures will be applied in a straightforward fashion to investigate i) the L5-dependence 
of the electron density function in a Hartree-Fock approximation; ii) the density 
functions of the atoms in their ground state, throughout the periodic table, based on 
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the density function alone; and iii) a quantization of relativistic effects by comparing 
density functions from Hartree-Fock and Dirac-Hartree-Fock models. 



7.1 On the LS-term dependence of atomic electron density 
functions 



From the developments on the density function in section 3.2 it is clear that the 
L5-dependent restricted Hartree-Fock approximation yields L5-dependent Hartree- 
Fock equations for atoms with open sub shells. In Hartree-Fock this dependence 
can be traced back to the term-dependency of the coulomb interaction. In the single 
incomplete shell case, corresponding to the ground state configurations we are in- 
terested in for the present study, the term dependency is usually fairly small. Froese- 
Fischer compared mean radii of the radial functions 1 17|. The differences in (r) for 
the outer orbitals between the values obtained from a Hartree-Fock calculation on 
the lowest term and those for the average energy of the configuration are of the or- 
der of 1-5 %. Although the L5-dependency does not show up explicitly in the direct 
and exchange potentials of the closed-subshell radial Hartree-Fock equations, the 
closed-subshell radial functions are ultimately L5-dependent through the coupling 
between the orbitals in the HF equations to be solved in the iterative procedure, but 
these relaxation effects turn out to be even smaller. 



As explicitly indicated through ( 3 1 1 the density built from the one-electron ra- 
dial functions could therefore be L5-dependent but this issue has not yet been inves- 
tigated quantitatively. Combining the term-dependent densities = Pu^LaSa ^^'^ 
Pb = PubLbSb for the same atom in the same electronic ground state configuration, 
but possibly different states (and adopting the Dirac 5-function for Q) for eval- 



uating the quantum similarity measure Zab of ( 68 1, the similarity matrix can be 
constructed according to (69 1. Its matrix elements have been estimated for the np^ 



configuration of Carbon (« = 2) and Silicon (« — 3) |6 |. As expected, the deviation 
of the off-diagonal elements from 1 is very small, the HF orbitals for the different 
terms ^P, 'D and '5 being highly similar, although not identical. 



7.2 A study of the periodic table 



As a first step to the recovery of the periodic patterns in Mendeleev's table, Carbo's 
quantum similarity index ( |69| was used, with the Dirac-5 as separation operator. In 
this case the expression (|69|l reduces to an expression for shape functions ([66]l. 



For the evaluation of the QSI in expression ( 69 1, we used atomic density func 



tions of atoms in their ground state e.g. corresponding to the lowest energy term. As 



elaborated in section 3.2 the involvement of all degenerate magnetic components 



allows to construct a spherical density function. For the density functions in this 
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Study, we limited ourselves to the Hartree-Fock approximation where no correla- 
tion effects are involved and the state functions are built with one CSF. 

In figure]?] we extract, as a case study from the complete atom QSI-matrix, the 
relevant information for the noble gases. Here the similarities were calculated using 
the Dirac delta function as separation operator From these data it is clear that the 
similarity indices are higher, the closer the atoms are in the periodic table (smallest 
AZ,Z being the atomic number). The tendency noticed by Robert and Carbo in 
ll45l is regained in the present study at a more sophisticated level. It can hence be 
concluded that the QSI involving p(r) and evaluated with 5(r — r') as separation 
operator Q., does not generate periodicity. 




Z 



Fig. 4 Quantum similarity indices for noble gases, using the Dirac-delta function as separation 
operator. 



The discussion of the work on the retrieval of periodicity is continued below 
in sections |9.1| and |9.2| where concepts from information theory are employed to 
construct a functional which quantifies the difference between two density functions 
in a different way. 



7.3 On the influence of relativistic effects 

In this section we investigate the importance of relativistic effects for the electron 
density functions of atoms. From the relativistic effects on total energies one can 
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infer these effects have implications for the electron densities. The effect of relativ- 
ity on atomic wave function has been studied in the pioneering work of Burke and 
Grant [46 1 who presented graphs and tables to show the order of magnitude of cor- 
rections to the hydrogenic charge distributions for Z ~ 80. The relative changes in 
the binding energies and expectation values of r due to relativistic effects are known 
from the comparison of the results obtained by solving both the Schrodinger and 
Dirac equations for the same Coulomb potential. The contraction of the n.s-orbitals 
is a well known example of these relativistic effects. But as pointed out by De- 
sclaux in his ''Tour historique" f47\, for a many-electron system, the self-consistent 
field effects change this simple picture quite significantly. Indeed, contrary to the 
single electron solution of the Dirac equation showing mainly the mass variation 
with velocity, a Dirac-Fock calculation includes the changes in the spatial charge 
distribution of the electrons induced by the self-consistent field. 

We first illustrate the difference of the radial density functions D{r) defined as 



(see also expression ( 34 1) 



D{r)=47Zr^p{r) , (70) 

calculated in the Hartree-Fock (HF) and Dirac-Fock (DF) approximations for the 
ground state 6p^ ^Pq of Pb I (Z 82) according to equations (65 1 and (70 1, respec- 



tively. These are plotted in figure |5j which shows the global relativistic contraction 
of the shell structure. 




oi < ' < ' < 1 ' ' < 1 

0.2 0.4 0.6 0.8 1 

r/an 

Fig. 5 DF and HF density distributions D{r) = 4;rr^p(r) for the neutral Pb atom (Z = 82). The 
contraction of the first shells is clearly visible. 
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Employing the framework of QSI to compare non-re lativistic Hartree-Fock elec- 
tron density functions p^^(r) with relativistic Dirac-Fock electron density func- 
tions p^^ (r) for a given atom, the influence of relativistic effects on the total density 
functions of atoms can be quantified via the QSI defined as 



ZHF,DF(5)=/p"F(r)5(r-r')pi^P(r')^/r^/r' 

Zhf. df(^) 



SIs = 



^/Zhf.df(5)^Z„f.bf(5) 



(71) 
(72) 



where 5 is the Dirac-5 operator. 

In figure|6]we supply the QSI between atomic densities obtained from numerical 
Hartree-Fock calculation and those obtained from numerical Dirac-Fock calcula- 
tions, for all atoms of the periodic table. 

The results show practically no relativistic effects on the electron densities for 
the first periods, the influence becoming comparatively large for heavy atoms. To 
illustrate the evolution through the table the numerical results of the carbon group 
elements are highlighted in the graph in figure|6] From the graph it is also noticeable 
that the relativistic effects rapidly gain importance for atoms heavier than Pb (Z = 
86). 



C 0.94 



0.92 



0.9 



100 



Fig. 6 Similarity of non-relativistic Hartree-Fock with relativistic Dirac-Fock atomic density 
functions with highUghted resuks for the C group atoms (C, Si, Ge, Sn, Pb). 
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8 Analyzing atomic densities: concepts from Information Theory 

Nowadays density functional theory (DFT) is the most widely used tool in quan- 
tum chemistry. Its relatively low computational cost and the attractive way in which 
chemical reactivity can be investigated made it a good alternative to traditional wave 
function based approaches. DFT is based on the Hohenberg-Kohn theorems 1 1]. In 
other words an atom's or a molecule's energy - and in fact any other physical prop- 
erty - can be determined by evaluating a density functional. However, the construc- 
tion of some functionals corresponding to certain physical property, has proven very 
difficult. Moreover, to the present day no general and systematic way for construct- 
ing such functionals has been established. Although energy functionals, which are 
accurate enough for numerous practical purposes, have been around for some time 
now, the complicated rationale and the everlasting search for even more accurate 
energy functionals are proof of the difficulties encountered when constructing such 
functionals. In the domain of conceptual DFT, where chemical reactivity is investi- 
gated, a scheme for the construction of functionals, based on derivatives of the en- 
ergy with respect to the number of electrons and/or the external potential, has proven 
very successful f40', l48j. Inspiration for the construction of chemically interesting 
functionals has also come from information theory and statistical mathematics. The 
functionals used for analyzing probability distributions have been successfully ap- 
plied to investigate electron density functions of atoms and molecules. In this chap- 
ter we introduce those functionals and discuss several studies where they have been 
applied to construct chemically interesting functionals. 

Shannon is generally recognized as one of the founding fathers of information 
theory. He defined a measure for the amount of information in a message and based 
on that, he developed a mathematical theory of communication. His theory of com- 
munication is concerned with the amount of information in a message rather than the 
information itself or the semantics. It is based on the idea that - from the physical 
point of view - the message itself is irrelevant, but its size is an objective quantity. 
Shannon saw his measure of information as a measure of uncertainty and referred 
to it as an entropy. Since Shannon's seminal publication in 1948 1,49 1, informa- 
tion theory became a very useful quantitative theory for dealing with problems of 
transmission of information and his ideas found many applications in a remarkable 
number of scientific fields. The fundamental character of information, as defined 
by Shannon, is strengthened by the work of Jaynes ll50ll5Tl . who showed that it is 
possible to develop a statistical mechanics on the basis of the principle of maximum 
entropy. 

In the literature the terms entropy and information are frequently interchanged. 
Arih Ben-Naim, the author of " Farewell to Entropy: Statistical Thermodynamics 
Based on Information " 152] insists on going one step further and motivates "not 
only to use the principle of maximum entropy in predicting the probability distri- 
bution [which is used in statistical physics], but to replace altogether the concept of 
entropy with the more suitable information". In his opinion "this would replace an 
essentially meaningless term [entropy] with an actual objective, interpretable phys- 
ical quantity [information]". We do not intend to participate in this discussion at 
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this time, since the present chapter is not concerned with the development of in- 
formation theory itself, but rather with an investigation of the applicability of some 
concepts, borrowed from information theory, in a quantum chemical context. The 
interested reader can find an abundance of treatments on information theory itself 
and its applications to statistical physics and thermodynamics in the literature. 

Information theoretical concepts found their way into chemistry during the sev- 
enties. They were introduced to investigate experimental and computed energy dis- 
tributions from molecular collision experiments. The purpose of the information 
theoretical approach was to measure the significance of theoretical models and con- 
versely to decide which parameters should be investigated to gain the best insight 
into the actual distribution. For an overview of this approach to molecular reaction 
dynamics, we refer to Levine's work f53\. Although the investigated energy dis- 
tributions have little relation with electronic wave functions and density functions, 
the same ideas and concepts found their way to quantum chemistry and the chemi- 
cal reactivity studies which are an important study field of it. Most probably this is 
stimulated by the fact that atoms and molecules can be described by their density 
function, which is ultimately a probability distribution. The first applications of in- 
formation theoretical concepts in quantum chemistry, can be found in the literature 
of the early eighties. The pioneering work of Sears, Parr and Dinur |54| quickly lead 
to more novel ideas and publications. Since then, many applications of information 
theoretical concepts to investigate wave functions and density functions, have been 
reported. In ||55]| Gadre gives a detailed review of the original ideas behind and 
the literature on "Information Theoretical Approaches to Quantum Chemistry". To 
motivate our work in this field we paraphrase the author's concluding sentence: 

"Thus it is felt that the information theoretical principles will continue to serve as powerful 
guidelines for predictive and interpretive purposes in quantum chemistry." 

The initial idea in our approach was to construct a density functional, which re- 
veals chemical and physical properties of atoms, since the periodicity of the Table 
is one of the most important and basic cornerstones of chemistry. Its recovery on 
the basis of the electron density alone can be considered a significant result. In an 
information theoretical context, the periodicity revealing functional can be inter- 
preted as a quantification of the amount of information in a given atom's density 
function, missing from the density function of the noble gas atom which precedes 
it in the periodic table. The results indicate that information theory offers a method 
for the construction of density functionals with chemical interest and based on this 
we continued along the same lines and investigated if more chemically interesting 
information functionals could be constructed. 

In the same spirit, the concept of complexity has been taken under consideration 
for the investigation of electron density functions. Complexity has appeared in many 
fields of scientific inquiry e.g. physics, statistics, biology, computer science and 
economics |56 1. At present there does not exist a general definition which quantifies 
complexity, however several attempts have been made. For us, one of these stands 
out due to its functional form and its link with information theory. 
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Throughout this chapter it becomes clear that different information and complex- 
ity measures can be used to distinguish electron density functions. Their evaluation 
and interpretation for atomic and molecular density functions gradually gives us 
a better understanding of how the density function carries physical and chemical 
information. This exploration of the density function using information measures 
teaches us to read this information. 

Before going into more details about our research several concepts should be 
formally introduced. For our research, which deals with applications of functional 
measures to atomic and molecular density functions, a brief discussion of these mea- 
sures should suffice. The theoretical sections are followed by an in depth discussion 
of our results. In the concluding section we formulate general remarks and some 
perspectives. 



8.1 Shannon's measure: an axiomatic definition 

In 1948 Shannon constructed his information measure - also referred to as "entropy" 
- for probability distributions according to a set of characterizing axioms |49|. A 
subsequent work showed that, to obtain the desired characterization. Shannon's 
original four axioms should be completed with a fifth one |57|. Different equiva- 
lent sets of axioms exist which yield Shannon's information measure. The original 
axioms, with the necessary fifth, can be found in ll58l . Here we state the set of ax- 
ioms described by Kinchin ^55] [59] . 

For a stochastic event with a set of n possible outcomes (called the event space) 
{Ai ,A2, . . . ,A„} where the associated probability distribution P = {Pi,P2, . . . ,Pn} 
with Pi > for all / and Y!i=i Pi = the measure S should satisfy: 

1 . the entropy functional 5 is a continuous functional of P; 

2. the entropy is maximal when P is the uniform distribution i.e. P, = 1 /n; 

3. the entropy of independent schemes are additive i.e. S{Pa +Pb) = S{Pa) +S{Pb) 
(a weaker condition for dependent schemes exists); 

4. adding any number of impossible events to the event space does not change the 
entropy i.e. , P2, • • ■ , i',,, 0, 0, . . . , 0) = , P2 , ■ • • , i^,,) ■ 



It can be proven 11591 that these axioms suffice to uniquely characterize Shannon's 
entropy functional 

S ^ -kJ^P.logPi , (73) 
with k a positive constant. The sum runs over the event space i.e. the entire prob- 



ability distribution. In physics, expression (73 1 also defines the entropy of a given 
macro-state, where the sum runs over all micro-states and where f, is the probability 
corresponding to the i-th micro-state. The uniform distribution possesses the largest 
entropy indicating that the measure can be considered as a measure of randomness 
or uncertainty, or alternatively, it indicates the presence of information. 
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When Shannon made the straightforward generalization for continuous probabil- 
ity distributions P{x) 

S[P{x)] = -kJ P{x) logP(x) dx , (74) 

he noticed that the obtained functional depends on the choice of the coordinates. 
This is easily demonstrated for an arbitrary coordinate transformation y ~ g{x), by 
employing the transformation rule for the probability distribution p{x) 

q{y)=p{x)J-' (75) 

and the integrandum 

dy = Jdx , (76) 

where J is the Jacobian of the coordinate transformation and J ' its inverse. The 
entropy hence becomes 

/ ?(3')log(^W) dy^ j p{x)\og{p{x)J-') dx , (77) 

where the residual J^^ inhibits the invariance of the entropy. Although Shannon's 
definition lacks invariance and although it is not always positive, it generally per- 
forms very well. Moreover, its fundamental character is emphasized by Jaynes's 
maximum entropy principle, which permits the construction of statistical physics, 
based on the concept of information llSOllSTI . In the last decade several investiga- 
tions of the Shannon entropy in a quantum chemical context have been reported. 
Those relevant to our research are discussed in more detail below. 



8.2 Kullback-Leibler missing information 

Kullback-Leibler's information deficiency was introduced in 1951 as a generaliza- 
tion of Shannon's information entropy [8J. For a continuous probability distribution 
P{x), relative to the reference distribution Po{x), it is given by 

AS[Pix)\Po{x)] = J P^log ^ dx . (78) 



As can easily be seen from expression ( 77 1, the introduction of a reference prob 



ability distribution Pq{x) yields a measure independent of the choice of the coor- 
dinate system. The Kullback-Leibler functional quantifies the amount of informa- 
tion which discriminates P{x) from Poix). In other words, it quantifies the distin- 
guishability of the two probability distributions. Sometimes it can be useful to see 
i45'[P(.5t)|/b(-^)] as the distance in information from Pq to P, although strictly speak- 
ing the lack of symmetry under exchange of P{x) and Po{x) makes it a directed 
divergence. 
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Kullback-Leibler's measure is an attractive quantity from a conceptual and for- 
mal point of view. It satisfies the important properties positivity, additivity, invari- 
ance, respectively: 

1. AS[P{x)\Po{x)]>0; 

2. AS[P{x,y)\Po{x,y)] = AS[P{x)\Pq{x)] + AS[P(j)\PQ{y)] for independent events 
i.e.P{x,y) =P{x)P{y); 

3. AS[P{y)my)]^ AS[P{x)\Po{x)] ify = /(x). 

Besides the lack of symmetry, the Kullback-Leibler functional has other formal 
limitations e.g. it is not bound, nor is it always well defined. In [60| the lack of these 
properties was addressed and the Jensen-Shannon divergence was introduced as a 
symmetrized version of Kullback-Leibler's functional. In 161] the Jensen-Shannon 
distribution was first proposed as a measure of distinguishabiUty of two quantum 
states. Chatzisavvas et al. investigated the quantity for atomic density functions |62|. 

For our investigation of atomic and molecular density functions, as carrier of 
physical and chemical information, we constructed functional based on the defi- 
nition of information measures. In sections |9T| below, the research is discussed in 
depth. 



9 Examples from information Theory 

9.1 Reading chemical information from the atomic density 
functions 

This section contains a detailed description of our research on the recovery of the 
periodicity of Mendeleev's Table. The novelty in this study is that we managed to 
generate the chemical periodicity of Mendeleev's table in a natural way, by con- 



structing and evaluating a density functional. As discussed before in section 7.2 the 
comparison of atomic density functions on the basis of a quantum similarity index 
(using the 5{r\ —V2) operator), masks the periodic patterns in Mendeleev's table. 
On the other hand, the importance of the periodicity, as one of the workhorses in 
chemistry, can hardly be underestimated. Due to the Hohenberg-Kohn theorems, 
the electron density can be considered as the basic carrier of information, although, 
for many properties it is unknown how to extract the relevant information from the 
density function. This prompted us to investigate whether the information measures, 
which gained a widespread attention by the quantum chemical community, could be 
used to help extract chemical information from atomic density functions in general 
and help to regain chemical periodicity in particular 



Tempted by the interpretation of the Kullback-Leibler expression (78 1 as a tool 
to distinguish two probability distributions, the possibility of using it to compare 
atomic density functions is explored. To make a physically motivated choice of the 
reference density Pq{x) we consider the construction of Sanderson's electronegativ- 
ity scale ll63l . which is based on the compactness of the electron cloud. Sanderson 
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introduced a hypothetical noble gas atom with an average density scaled by the num- 
ber of electrons. This gives us the argument to use renormalized noble gas densities 



as reference in expression (78 i. This gives us the quantity 



ASP, = AS[pA{r)\poir)] = / Pa (r) log ^ dr , (79) 



. gA(r) 
(^o(r) 

where (r) and ox (r) are the density and shape function of the investigated system 
and C7o(r) the shape function of the noble gas atom preceding atom A in Mendeleev's 
table. The evaluation of this expression for atoms He through Xe shows a clear 
periodic pattern, as can be seen in Figure |7] 

Reducing the above expression to one that is based on shape functions only, leads 

to 

AS^ = AS[aA{r)\ao{r)] = ^ oa (r) log t/r (80) 

and its evolution is shown in Figure [8] The periodicity is clearly present and this 
with the fact that the distance between points in a given period is decreasing grad- 
ually from first to fourth row is in agreement with the evolution of many chemical 
properties throughout the periodic table. One hereby regains one of the basic charac- 
teristics of the Periodic Table namely that the evolution in (many) properties through 
a given period slows down when going down in the Table. The decrease in slope of 
the four curves is a further illustration. 



CO 



Fig. 7 Kullback-Leibler information p9\ versus Z for atomic densities with the noble gas of the 
previous row as reference. 
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Fig. 8 KuUback-Leibler infomiation (equation |80j) versus Z for atomic shape functions with the 
noble gas of the previous row as reference. 



9.2 Information theoretical QSI 

Continuing the search for periodic patterns based on similarity measures, as intro- 
duced in section |7] and motivated by the resuhs obtained in an information theoreti- 
cal framework in section [9T1 we will now combine the ideas from both efforts and 
construct an information theoretical similarity measure. 

For the construction of the functional in the above section, the choice to set the 
reference (the prior) density to that of a hypothetical noble gas atom, in analogy to 
Sanderson's electronegativity scale, was motivated and the particular choice lead to 
results which could be interpreted chemically. Following these findings one can see 
that it would be interesting to compare the information entropy, evaluated locally as 

ASP ir)^p^ir)\og^^^ (81) 

for two atoms by use of a QSM, which can be constructed straightforwardly, by 
considering the overlap integral (with Dirac-5 as separation operator) of the local 
information entropies of two atoms A and B 

ZAB{d)= f pA(r)log^|^p5(r)log^|^«'r. (82) 

A QSI can be defined by normalizing the QSM as before, via expression ( [69] l. The 
QSM and the normalized QSI give a quantitative way of studying the resemblance 
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in the information carried by the valence electrons of two atoms. The obtained QSI 
trivially simplifies to a shape based expression 



51(5) - 



J AS^{r)AS^{r)dr 

A S- {r)AS- {v)dr^jAS% {r)AS% {r)dr 



(83) 



To illustrate the results we select the QSI ( |83] l with the top atoms of each column 
as prior Formulated in terms of KuUback-Leibler information discrimination the 
following is evaluated. For instance, when we want to investigate the distance of the 
atoms Al, Si, S and CI from the N-column (group Va), we consider the information 



theory based QSI in expression (83 i, where the reference densities po and pg/ are 
set to Pn, Pa to Pai, psi, pp, etc. respectively and Pb to p/>, i.e. we compare the 
information contained in the shape function of N to determine that of P, with its 
information on the shape function of Al, Si, S, CI. Due to the construction a 1. is 
yielded for the element P and the other values for the elements to the left and to the 
right of the N-column decrease, as shown in figure |9] This pattern is followed for 
the periods 3 up to 6, taking As, Sb and Bi as reference, with decreasing difference 
along a given period in accordance with the results above. Note that the difference 
from 1 . remains small, due to the effect of the renormalization used to obtain the 
QSI. 



0.996 



0.994 



CO 




Fig. 9 Results of the information theory based QSI with the atom on top of the column as prior. 
The symbol in the legend indicates the period of the investigated atom and the nuclear charge 
Z-axis indicates the column of the investigated atom. (For example Ga can be found as a square 
Z = 5). 
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Results on the investigation of atomic density functions are reviewed. First, ways 
for calculating the density of atoms in a well-defined state are discussed, with par- 
ticular attention for the spherical symmetry. It follows that the density function of an 
arbitrary open shell atom is not a priori spherically symmetric. A workable defini- 
tion for density functions within the multi-configuration Hartree-Fock framework is 
established. By evaluating the obtained definition, particular influences on the den- 
sity function are illustrated. A brief overview of the calculation of density functions 
within the relativistic Dirac-Hartree-Fock scheme is given as well. 

After discussing the definition of atomic density functions, quantum similarity 
measures are introduced and three case studies illustrate that specific influences on 
the density function of electron correlation and relativity can be quantified in this 
way. Although no periodic patterns were found in Mendeleev's table, the methodol- 
ogy is particularly successful for quantifying the influence of relativistic effects on 
the density function. 

In the final part the application of concepts from information theory is re- 
viewed. After covering the necessary theoretical background a particular form of the 
Kullback-Liebler information measure is adopted and employed to define a func- 
tional for the investigation of density functions throughout Mendeleev's Table. The 
evaluation of the constructed functional reveals clear periodic patterns, which are 
even further improved when the shape function is employed instead of the density 
functions. These results clearly demonstrate that it is possible to retrieve chemically 
interesting information from the density function. Moreover the results indicate that 
the shape function further simplifies the density function without loosing essential 
information. The latter point of view is extensively treated in 1641 . where the authors 
elaborately discuss "information carriers" such as the wave function, the reduced 
density matrix, the electron density function and the shape function. 
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